Global implicit solver for multiphase multicomponent flow in porous media with multiple gas components and general reactions

Abstract In order to study the efficiency of various forms trapping including mineral scenarios for CO 2 storage behavior in deep layers porous media, highly nonlinear coupled diffusion-advection-reaction partial differential equations (PDEs) kinetic and equilibrium reactions modeling miscible multiphase multicomponent flow have be solved. We apply globally fully implicit PDE reduction method (PRM) developed 2007 by Kräutle Knabner one-phase flow, which was extended 2019 case two-phase with a pure gas Brunner Knabner. extend an arbitrary number gases gaseous phase, because is not only that threats climate, usually accompanied other climate killing gases. The application PRM leads equation system consisting PDEs, ordinary equations, algebraic equations. Finite Element discretized / Volume stabilized are separated into local global but nevertheless resolution function evaluated aid nested Newton solver, so our solver implicit. For phase disappearance, we use persistent variables lead semismooth formulation solved method. present injection mixture layers, investigate change effects context gases, technique. technical framework also applies fields such as nuclear waste or oil recovery.